Survival probability of a particle in a sea of mobile traps: A tale of tails 
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O 

We study the long-time tails of the survival probability P(t) of an A particle diffusing in d-dimensional 
£SJ ' media in the presence of a concentration p of traps B that move sub-diffusively, such that the mean square 

displacement of each trap grows as t 7 with < 7 < 1. Starting from a continuous time random walk (CTRW) 
description of the motion of the particle and of the traps, we derive lower and upper bounds for P(t) and show 
that for 7 < 2/(d + 2) these bounds coincide asymptotically, thus determining asymptotically exact results. The 
asymptotic decay law in this regime is exactly that obtained for immobile traps. This means that for sufficiently 
subdiffusive traps, the moving A particle sees the traps as essentially immobile, and Lifshitz or trapping tails 
remain unchanged. For 7 > 2/(d + 2) and d < 2 the upper and lower bounds again coincide, leading to a decay 
law equal to that of a stationary particle. Thus, in this regime the moving traps see the particle as essentially 
immobile. For d > 2, however, the upper and lower bounds in this 7 regime no longer coincide and the decay 
■ law for the survival probability of the A particle remains ambiguous. 

o , 

PACS numbers: 02.50.-r,82.40.-g,89.75.Da 



I. INTRODUCTION 

> 

In the early 1960's, I.M. Lifshitz published his famous analysis of the low-energy tails of the density N(E) of states of an 
electron in a medium with randomly scattered immobile impurities lUlH]. He demonstrated that, in <i-dimensions, the spectrum 
Q N(E) of the random Schrodinger operator decays as 



N(E) ~ exp (-const E- d/2 \ , E -> 0. (1) 



a 

o 

This exponential decay is in striking contrast to the polynomial decay associated with a periodic Schrodinger operator, and is 

caused by the presence of arbitrarily large, albeit rare, regions without impurities. 
, A decade later Balagurov and Vaks [3] and Donsker and Varadhan [4] published their celebrated work on trapping kinetics 
■ in a medium with randomly placed immobile traps (called B). They found that the long-time decay of the survival probability 
' P B (t) of a particle (called A) diffusing in such a d dimensional medium follows the stretched-exponential law 

(N 

i P B (t) ~ exp (-k dP ^^\D A t) d l^) , (2) 

where p is the mean density of traps, Da is the particle diffusion coefficient, and kd is the d-dependent constant (01 



d + 2 



k d =[^)(dv d f^(2 Z jf {d+2) . (3) 



Here Zd is the first zero of the Bessel function J(d-2)/2{ z ) an d v d = 2n d ^ 2 /[dT(d/2)] denotes the volume of a d-dimensional 
sphere of unit radius. The superscript B on the survival probability emphasizes that the traps B are immobile. Balagurov and 
Vaks obtained the decay form in Eq. (|2]) exactly for d = 1. While they furthermore deduced this behavior for general d by 
noticing the close mathematical connection between the trapping problem and the Lifshitz problem, Donsker and Varadhan were 
the first to obtain the rigorous exact solution of this essentially many-body problem in d dimensions. They determined not 
only the time dependence but also the decay coefficient kd of Eq. ©. Note that this law stands in stark contrast to the purely 
exponential decay predicted by standard chemical kinetics for the reaction A + B — > B of mobile particles and traps, and even 
to the predictions of the Smoluchowski approach based on a reaction-diffusion equation, 

{exp[- j o(4L>t/ 7 r) 1 / 2 ], d=l 
exp[~4TrpDt/\n(Dt/a 2 )}, d = 2 (4) 
exp(-4TT P Dt), d = 3 

where D is the sum of the diffusion coefficients of the particle and the traps |Ht]. 
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In mathematical analogy with the source of the Lifshitz tail, the decay law (O arises from the presence of arbitrarily large 
regions without traps in which the particle can diffuse for a relatively long time before being trapped. Although the survival 
probability of a particle in any one such region is purely exponential, the ensemble average over a random distribution of such 
trap-free voids produces the anomalously slow decay. 

The decay law in Eq. (f2]l has been also generalized to trapping reactions on fractals and in inhomogeneous structures @], 
to A particles performing subdiffusive motion 0, [R & O or attached to the extremities of polymer chains ifTTIl . and also to 
agglomerated random distributions of traps distributed on immobile polymer chains in solution lfl2ll or in clusters Q. The 
history of this problem and many other results have been summarized in several reviews (see, e.g.. lll4l,[l5ll ). 

The survival probability Eq. (O [as well as the Lifshitz tail result of Eq. (fl]i] is valid only when the traps (or impurities) are 
strictly immobile. Indeed, if one allows them to diffuse, no matter how small the diffusion coefficient, the particle survival 
probability is described by a faster decay law. As proved by Bramson and Lebowitz lllql (see also \ Y£\), when both species 
diffuse, the survival probability of the A particle at long times is given by 

!exp (-pcit 1 / 2 ) , d = 1 
exp(-pc 2 t/]nt), d=2 (5) 
exp (— po$t) , d = 3. 

The time dependences are the same as in the Smoluchowski problem, but the constants ci , C2 and C3 are in general different, 
may depend on the diffusion coefficients of the particle and of the traps (but see below), and were not determined in this original 
work. 

When traps diffuse, the many-body trapping effects captured by Eq. © are thus no longer applicable and the A particle 
survival probability decays with time according to the faster time dependencies in Eq. ©. The underlying fluctuation mechanism 
governing the trapping dynamics has changed when the traps are allowed to move. While the time dependence in the decay 
laws © is consistent with Smoluchowski-like results, which in fact represent a two-body approximation, the decay amplitudes 
Cd are not simply functions of the sum of the diffusion coefficients. In particular, the exact form (including the coefficients) of 
the leading large-i behavior for d — 1 and d = 2 has only recently been found by Bray and Blythe [18] and, surprisingly, ci 
and C2 depend only on the diffusion coefficient Db of the traps and are independent of Da, the diffusion coefficient of the A 
particle! This implies that the survival probability of the diffusing A particle in this scenario is asymptotically identical to that 
of an A particle that remains still, that is, to P A {t). 

This remarkable result has subsequently been extended to systems in which both the particles and the traps move subdiffu- 
sively lp7L fsL f^L fTotl . In particular, it is again found that the survival probability of a subdiffusive A particle in a sea of subdiffusive 
traps is identical to that of and A particle that remains still (as in the Bramson-Lebowitz scenario), P(t) ~ P A (t). HoweverJt is 
still not clear what happens when traps move subdiffusively and the particle is diffusive. In one dimension it was shown in Qlal 
that the asymptotic survival probability of the diffusing particle is the same when surrounded by subdiffusive traps characterized 
by a mean square displacement that grows as V with 7 < 1 as it is for a stationary particle, provided the subdiffusive traps move 
sufficiently rapildy (7 > 2/ 3), again akin to the Bramson-Lebowitz scenario. However, few results seem to be currently available 
for d > 1 and/or if the subdiffusive traps are "too slow." An interesting question is then the following: will the trapping tails 
associated with immobile traps withstand sufficiently slow sub-diffusion, thus leading to a survival probability such as that of 
the Donsker-Varadhan result, or will they switch to appropriately rescaled Bramson-Lebowitz forms? This is the main question 
that we pose in this paper. 

In this pursuit, we follow the general idea HHl of constructing lower and upper bounds on the survival probability of a diffusive 
particle in a sea of subdiffusive traps whose motion is described by a continuous time random walk that starts at time t — lfl9ll . 
If these bounds coincide asymptotically, then we can extract an exact asymptotic result for the survival probability of the particle. 
We write the asymptotic survival probability in the form 

fexp(-r), d^2 
1 ' ~ 1 exp (-0*7 In t) d=2, 

where the constants 8 and z depend on 7 and on dimension, and explore whether a convergence of upper and lower bounds 
provides the exponents z and perhaps even the exponential prefactor 9. We broadly anticipate our results by noting that if the 
subdiffusive particles are sufficiently slow, specifically if 7 < 2/ (d + 2), these bounds lead to a Donsker-Varadhan behavior, and 
the Lifshitz or trapping tails thus remain unchanged. If 7 > 2 /(d+2) and d < 2, then the bounds lead to a Bramson-Lebowitz 
behavior that is, a survival probability behavior associated with a stationary particle (note that, albeit for a different system, this 
result is implicit in 1 1 Offl ; that work addresses diffusive particles in a fractal medium, while here we are considering subdiffusive 
motion in Euclidean geometries). However, if 7 > 2/(d + 2) and d > 2 we are not able to establish a unique asymptotic 
behavior. The detailed results are exhibited later in the paper. 

In Sec. we formulate the model of moving particle and traps. In Sec. [Ill] we calculate a lower bound on the survival 
probability of the particle, and in Sec. [IV] we obtain two upper bounds. The consequences of these bounds on the survival 
probability are collected in Sec. [V] A brief recapitulation of the results is given in the concluding section. 
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II. THE MODEL 



In this section we formulate the general model that allows us to highlight the approximations made to obtain upper and 
lower bounds on the survival probability. Consider a d-dimensional system of volume V containing a single diffusive A particle 
of radius a and K randomly moving point traps B (a finite trap radius would add nothing interesting to the problem for our 

purposes). The initial position of A is the origin, and one realization of its trajectory is denoted by a t . The starting points 

(k) 

-Bq , k = 1, 2, • • • K of the traps are randomly (Poisson) distributed throughout the volume, and a trajectory of the fcth trap is 
denoted by B { t k) = B { Q k) + b { t k) . 
Next we define 

J oo, \x\ < a, 

v(xj = < (7) 
1 0, otherwise. 

The indicator function ^ [a t , {b[ k ^ }] of the event that the A particle has not met any of the Bs up to time t for a given realization 
of their trajectories can then be written as 



K K 

-b^-B^dr 



nau{b? ) }] = \{naub ( t ) ] = TT exp "/ »K " 4' 

k =i k=i L Jo 



(8) 



Consequently the A particle survival probability is 

P(t)=E^Sj[ (I jfdB«25« {*[««,&« ]})}, (9) 

where the symbol E^ {...} denotes an average over all A particle trajectories a t that start at the origin. The symbol E^ {...} 
denotes an average over all B trajectories b t whose starting point is Bq. We have labeled each of the latter trajectories by a 
trap label k, but since the traps move independently we can omit the label. Furthermore, we go to the thermodynamic limit 
K, V — > oo while keeping the ratio p = K/V fixed. This leads to the A survival probability 

P(t) = E^{e- f > E ' b) ( w ^- b ^}, (10) 

where W[a t — b t ] is the functional of the trajectories a t and b t 

(11) 



W[a t -b t ] = J dB (l - e -JS*(.°r-*T-Bo)*rj . 



The exact problem has thus been reduced to an effective two-body problem involving a single A particle and a single B trap. 
Nevertheless, it unfortunately does not seem possible to evaluate the survival probability exactly from this expression, the main 
mathematical difficulty being that the average over all possible trajectories b t has to be performed for each fixed a t . Only after 
this average has been performed can one then go on to carry out the further average over the A particle trajectories. This appears 
to be a non-tractable mathematical problem, and recourse has to be made to controllable approximations. We do it here by 
constructing lower and upper bounds on P(t) and identifying conditions and regimes where these converge asymptotically. 



III. LOWER BOUND 

A lower bound was originally devised in lfl8ll for diffusive particles and traps and was extended in J3, H, [Toll to the case of 
(sub)diffusive particles and traps. We adapt this method to the current situation in which the traps B perform continuous-time 
random walks, so that we may restate this bound in the language introduced above. Following Bray and Blythe, we introduce 
a notional spherical volume Vi of radius I centered at the origin, and pick only those realizations of initial trap distributions for 
which this volume is completely devoid of traps. The probability that the region will initially be empty of traps is 

P void {l) ~ exp (-v dP l d ) . (12) 

We furthermore introduce the probability t) that the A particle does not leave the notional volume V/ during time t lfl8tl . 



P A (l-t) - exp (-zjD A t/(l - a) 2 ) 



(13) 
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Thirdly, we introduce the probability that no B particle enters the notional volume up to time t, that is, the probability that an 
immobile rf-dimensional target of radius I survives up to time t in the presence of a concentration p of traps all of which perform 
subdiffusive motion |9[] : 



P B (l;t) 



{^D B t 1 ) d/2 

exp I ~ /, r(i-d/2)r(i + 7 d/2) 



exp ( -P^rrS^r I > d = 2 d4) 



r(l + 7)ln(4D B *70 
(_ 2{d~1)TT d / 2 l d - 2 D B t' 1 ' 

exp l p r(rf/2)r(i + 7 ) 




Here T is the gamma function, and D B is the anomalous diffusion coefficient of the traps, that is, the coefficient in the mean 
square displacement relation (r 2 ) = 2dD B t 1 /T(\ +7). 

Since the functional in Eq.© is positive definite, the two latter constraints on the trajectories of A and the Bs naturally lead 
to a lower bound on the survival probability P(t). Furthermore, for these constrained trajectories the functional in Eq.dTTb is 
strictly equal to zero, and hence Eq. (0 restricted in this way is equal to unity. As a consequence, the probabilities associated 
with the random processes a t and b t subject to these constraints can simply be factored, immediately leading to the lower bound 
onP(f) 

P(t) > P votd (l)P A (l;t)P B (l;t) = P L (l-t). (15) 

Finally, we note that this lower bound Pl (I ; t ) is in fact a family of lower bounds dependent on the radius I of the notional volume 
separating particle and traps. This radius can be chosen to give the best lower bound, that is, the maximal lower bound, which 
we simply denote as Pi(t). The optimal radii are shown explicitly in Appendix lAl and are shown to depend on dimensionality 
and on 7. The associated optimal lower bounds for the survival probability of A, which also depend on dimensionality and on 
7, then immediately follow. We thus have 

P{t) > P L (t) (16) 

where, for d < 2, 




2 

d+2 



-k(d)p 2 ^ d+2 \D A t) d ^ d +^) , 7 < 

- <! m ( pj^DB^fl 2 \ (17) 

r(l -d/2)r(l + 7 d/2) J ' d+2 ' 



For d — 2 we have 



'exp(-(4 7 rz 2 2 /0 ^t) 1 /2) ; 7 <| 



For d > 2 the results are 



eX P'- (27-l)r(l+7)lnt )' ^ 



2 

d+2 



( exp (-k{d)p 2 ^ d +V {D A t) d ^ d +V) , 7 < 

P ^ - c::p ( J {d-2), d ' 2 pD B t^ /d (4DAt\ d - 2)/d \ ^ (19) 

j^expi ay r ( rf/2 )r(l + 7 ) ) ) J> 7 > d+2- 

There is thus a change in behavior of the lower bound in all dimensions when 7 crosses the value 2/(d + 2). 



IV. UPPER BOUNDS 
A. Pascal Principle 

An upper bound on the survival probability P(t) of the diffusive particle was recently derived on the basis of the so-called 
"Pascal Principle," which states that the A particle survives longer if it remains still than if it moves. The problem with a 
static A particle and moving traps is the so-called "target problem," and consequently we label this upper bound on the survival 
probability as Py.tn.meAf)- The Pascal Principle was conjectured in IU8I1 and was proved in B20I1 for d < 2 for conventional 
diffusive motion. In Il2lll the statement was proved for a rather general class of random walks on d-dimensional lattices. A 
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similar statement was introduced more than a decade earlier in l22ll in the context of excitation energy migration. This upper 
bound is given by Eq. ( TT4l > obtained in [9] if we set I = a, the radius of the A particle. For visual ease we explicitly rewrite 
Eq. (TBI) with this replacement, 



P 



U, target 




{A-KD B f) d/2 

p r(i-d/2)r(i + 7 d/2) 

r(l + 7 ) In {lD B V/a i ) 
27r d ^a d - 2 D B t^ 



r(d/2)r(i + 7 ) 

To clearly tie together the various notations introduced so far, we note that 

P B (a;t) = P A (t) = P, 
and 

P(t) < Pu,target(t) 



U,target 



(0 



d < 2 
d= 2 
d > 2. 



(20) 



(21) 



(22) 



However, this upper bound when associated with the lower bound is not always sufficiently tight to provide the desired 
information about the asymptotic survival probability of the diffusive particle. We thus introduce and alternative new upper 
bound, which in some cases is lower than the above. In the next section we then explicitly pick the bounds to be used and exhibit 
the information that can be obtained from them. 



B. Anti-Pascal Principle 

This new upper bound is based on what we will call the "Anti-Pascal Principle." We will show that the worst possible strategy 
for traps in their search for a target is to remain immobile. Random motion, even uncorrected with the motion of the target, 
enhances the probability to encounter the target. In other words, the diffusing particle A survives longer if the traps remain still. 
The problem with a moving A particle in a sea of static traps is the so-called "trapping problem," and so we label this upper bound 
on the survival probability as Pu,tra PP ing- Again, to clearly tie together various notations we note that P B (t) = Pu,tra PP ing(t), 
and the upper bound just introduced then says that 

P(t) < Pu,tra P ping(t) = P B (t) = fiW {e^M} . (23) 

To prove Eq. (l23l we make use of Jensen's inequality for convex functions, which for our model can be stated as 

e -^ W {^-*.]}^^ 6) | e _^ [at _ bt] J (24) 

(Note that this inequality is generally used to derive a lower bound. Indeed, when applied to the average in Eq. ( fT0T > it 
yields a lower bound which is exactly the Smoluchowski-type result [17].) Consequently, we have the following upper bound: 

P(t) < {E^ { e -f> w l a *- b AX = E^ L-" W M\ , (25) 

where E 1 -^ {...} denotes an average with respect to the trajectories c t — a t ■ — b t of a "fictitious" particle C of radius a which 
starts its motion at B . 

Note that W[ct] has a clear geometric interpretaion - it defines the volume swept by the fictitious particle C within the time 
interval (0, t) and is thus an analog of the so-called Wiener sausage for conventional diffusive motion. Clearly, W[ct] is a non 
decreasing function of time t or, more precisely, of the number N of jumps made by the fictitious particle within the time interval 
(0, t). This number has two contributions, N = Na + N B , where Na is the number of jumps made by particle A and N B the 
number of jumps made by a trap. Setting N B = 0, i.e., supposing that the trap is immobile, clearly diminishes the total number 
of jumps and hence diminishes W[ct]. In this way we tighten the bound in Eq.d25il to arrive at the desired inequality ( 1231 We 
have thus established the second upper bound ( l23l for the survival probability of A. The maximal information about the survival 
probability is thus obtained from the lower bounds presented in Eqs. dT6b-(fT9T> and the smaller of the upper bounds d22b with d20b 
and ([23) with ©. 
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V. COLLECTING RESULTS AND BOUNDING THE SURVIVAL PROBABILITY 

In this section we collect our detailed results, exhibited by displaying values for the constants z and 9 in Eq. © when the 
upper and lower bounds converge asymptotically, and providing bounds for the exponent z when they do not. The results for the 
exponent z as a function of the subdiffusive trap exponent 7 are shown for integer dimensions d = 1,2, and 3 in Fig.Q] Even 
before displaying the explicit results we note the following: 

• When 7 < 2/{d + 2) the asymptotic survival probability of the diffusive A particle in the sea of subdiffusive traps is the 
same as it is for immobile traps. The question asked in the introduction, whether trapping tails associated with immobile 
traps withstand sufficiently slow sub-diffusion, is thus answered in the affirmative, with a precise dimension-dependent 
characterization of what is "sufficiently slow." This behavior, obtained by the asymptotic convergence of the bounds Pl (i) 
and Pjj,trapping (t), is shown for the exponent z by the thick solid lines with superimposed circles in Fig.Q] In this regime 
the survival probability is thus of the Donsker-Varadhan form. 

• When 7 crosses the value 2/(d + 2) there is a kind of dynamical phase transition. For d < 2 the survival probability of 
the diffusive particle in the sea of now more rapidly moving subdiffusive traps decays as it would if the particle remains 
immobile. This is indicated by the thin solid lines with superimposed circles in the first two panels of Fig. Q] and results 
from the asymptotic convergence of the bounds Pi(t) and Pu,target(t)- In this regime the survival probability is thus of 
the Bramson-Lebowitz form generalized to subdiffusive traps. 

• For d > 2 and 7 > 2/ (d+ 2) the situation is left somewhat uncertain: we are only able to bound the decay exponent but not 
determine it uniquely, because neither upper bound converges asymptotically to the lower bound. All we can say is that the 
survival probability decay exponent z lies in the triangular region bounded by the thin solid line with superimposed circles 
(associated with the lower bound), the thick solid line (associated with Pu,trapping), and the thin solid line (associated 
with Pu.target), as indicated in the figure caption. 

Finally, we collect the explicit results described above in Table [V] The reported results are for the exponents z and the 
prefactors 8 in Eq. (O. Again, we note that the results for z are sketched explicitly for d — 1, 2, 3 in Fig.Q] 

VI. CONCLUSIONS 

In conclusion, we have rigorously determined the survival probability of a particle diffusing in a d-dimensional medium in 
the presence of a concentration of traps performing subdiffusive random motion. We have arrived at our results by calculating 
a lower bound and two alternative upper bounds for the survival probability. One of the upper bounds, based on the "Pascal 
Principle," is obtained by assuming the particle to remain immobile. The other, based on the "Anti-Pascal Principle," is found 
by assuming that the traps remain immobile. We can then choose the tighter (lower) upper bound. Results for the asymptotic 
survival probability of the particle can thus be extracted if the lower bound and one of the upper bounds converge asymptotically. 




0.0 0.2 0.4 0.6 0.8 1.0 ' 0.00 0.25 0.50 0.75 1.00 0.0 0.2 0.4 0.6 0.8 1.0 

7 7 7 



FIG. 1: Asymptotic results for the survival probability of a diffusive particle in a sea of subdiffusive traps. Plotted is the exponent z as 
expressed in Eq. {6), or bounds on this exponent, as a function of the trap subdiffusion exponent 7. First panel: d = 1; second panel: d = 2. 
Third panel: d — 3. The solid lines are upper bounds; specifically, the thick solid lines are Pascal upper bounds (immobile particle) while the 
thin solid lines are Anti-Pascal upper bounds (immobile traps). The solid lines with superimposed circles are lower bounds. Thus the first and 
second panels show that for d < 2 the Pascal (Donsker-Varadhan) and lower bound lines coincide for 7 < 2/(d + 2). The behavior changes 
when 7 crosses the value 2/(d + 2), leading to a coincidence of the Anti-Pascal (Bramson-Lebowitz) and lower bound lines. For d < 2 the 
asymptotic exponent is thus obtained for all values of 7. The third panel for d = 3 shows convergence of the upper (Pascal) and lower bounds 
for 7 < 2/(d + 2) . However, when 7 > 2/(d + 2) = 2/5 the situation is uncertain. The exponent z now lies in the triangular region delimited 
by the lower bound and the Pascal (2/{d + 2) < 7 < d/(d + 2)) or Anti-Pascal (d/(d + 2) < 7 < 1) upper bound. 
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Dimension d 


Trap Subdiffusive 
Exponent 7 


Optimal 
Upper Bound 


Survival Probability 
Exponent z 


Survival Probability 
Prefactor 9 


1 < d < 2 


< 7 < -r-^ 

— ' — a + 2 


P\J ^trapping 
PlJ, tar get 


d 

d-\- 2 


k(d)p 2 ^ d+ ^D d ^ d+2) 

P (A7YD B ) d/2 


ih^^ 1 


r(i-d)r(i + 7 d/2) 


a — 2 


0<7< 2 
^<7<1 


, trapping 
PlJ, tar get 


1 

1 


k(2){pDA) ' 

1 ^ r(i + 7 ) 2 

7 - 4tv P D b - ^7 - 1 


d>2 


°^^7JT2 

2 ^ ^ d 
1+2^^1^1+2 

1+h^^ 1 


PlJ, trapping 
PlJ, tar get 
PlI, tar get 


d 

d + 2 

d ^ , d-2 , 27 

ir^^^^hr + it 
7<*<^ + ^ 


fc(d)p 2 /( d + 2 )/3l /(d+2) 
7 
? 



TABLE I: Collected results for the asymptotic survival probability exponent 2 and prefactor 9. 



Following this procedure, we have shown that when the dynamical exponent 7 characterizing the growth of the second moment 
of the displacement of the traps is less than 2/(d + 2), that is, if the traps move sufficiently slowly, the decay of the survival 
probability of the diffusing particle in any dimension is given exactly by the Donsker-Varadhan result obtained for immobile 
traps. When the traps move more rapidly than this, i.e., when 7 > 2/(d+ 2), then in dimensions d < 2 the survival probability of 
the particle is identical to that of a stationary A particle in a sea of mobile traps. For higher dimensions, d > 2, our results do not 
uniquely determine the survival probability of the particle in this 7 regime, but they do provide tighter bounds than previously 
known. We close by noting the well-known difficulties that may stand in the way of reaching and verifying some asymptotic 
results via numerical simulations (see [8] and references therein). Such simulations continue to pose an interesting challenge. 
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APPENDIX A: OPTIMAL NOTIONAL VOLUME 



In this appendix we explicitly show the optimal radius / that leads to the maximal lower bound on the survival probability of 
the diffusive particle in the presence of a sea of subdiffusive point traps. Combining the expressions in Eqs. (fT2l . ( fT3l and (fl4] > 
and differentiating their product with respect to I, we find that the optimal / depends on dimensionality and on 7. For d < 2 we 
find 



2zlD A t \ 1/(d+2) 
dv d p ) 



(Al) 
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For d = 2we obtain 



For d > 2 we find 



(zlD A t\ l/i 



1<\ 



l~ )\ n P / (A2) 
11 ' z 2 2 r(l + 1 )D A t \ 1/2 ( lGKpD%er \ 




9 j n ,\V(<i+2) 

2z d D At ) < _2_ 

^r(rf/2)r(i + 7 )^ 1 -^ 1/d _ 2 



(d-2) 2 n d ' 2 pD P / ' 7>d + 2 ' 
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